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ABSTRACT 

Only by incorporating various forms of feedback can theories of galaxy forma- 
tion reproduce the present-day luminosity function of galaxies. It has also been ar- 
gued that such feedback processes might explain the counter-intuitive behaviour of 

'downsizing' witnessed since rcdshifts z ~l-2. To examine this question, observations 
spanning 0.4 < z < 1.4 from the DEEP2/Palomar survey arc compared with a suite 
of equivalent mock observations derived from the Millennium Simulation, populated 
with galaxies using the Galform code. Although the model successfully reproduces 
the observed total mass function and the general trend of 'downsizing', it fails to ac- 
curately reproduce the colour distribution and type-dependent mass functions at all 
redshifts probed. This failure is shared by other semi-analytical models which collec- 
tively appear to "over-quench" star formation in intermediate-mass systems. These 
mock lightcones are also a valuable tool for investigating the reliability of the obser- 
vational results in terms of cosmic variance. Using variance estimates derived from 
the lightcones we confirm the significance of the decline since z ~ 1 in the observed 
number density of massive blue galaxies which, we argue, provides the bulk of the 
associated growth in the red sequence. We also assess the limitations arising from 
cosmic variance in terms of our ability to observe mass-dependent growth since z ^ 1. 



1 INTRODUCTION 

The physical picture of how galaxies assemble has changed 
markedly over the past decade. A pure 'hierarchical dark 
matter model' in which gas cooling and subsequent star for- 
mation occurs in synchronisation with the growth, via grav- 
itational instabihty, of their parent dark matter halos fails 
to reproduce the local luminosity function of galaxies (Ben- 
son et al. 2003; Somerville & Primack 1999; Kauffman et 
al. 1999) and has been challenged by the presence of mas- 
sive (~ 10^^ Mq) galaxies at redshifts z ~2 (Glazebrook 
et al. 2004; Cimatti et al. 2004; van Dokkum et al. 2006). 
As a result, a new paradigm has emerged which argues for 
the importance of 'feedback' processes that serve to govern 
the star formation rate in a galaxy. As the efficacy of these 
processes depends on the mass of the host galaxy, so it is 
possible to reconcile the predictions of the standard CDM 
model with the local galaxy luminosity function (Croton et 
al. 2006; Bower et al. 2006). 

Despite this progress, the physical basis of the feedback 
processes incorporated into the recent semi-analytic models 
remains largely untested. The most effective way of sup- 
pressing star formation and hence inhibiting further growth 



in massive galaxies is "radio mode" feedback (Croton et al. 
2006; Bower et al. 2006), where additional gas cooling in 
halos in which the cooling time is longer than the dynami- 
cal time is prevented by low levels of accretion onto central 
supermassive black holes. Bower et al. (2006) have argued 
that such a process can lead naturally to a characteristic 
mass scale associated with the transition between cooling 
on a hydrostatic timcscalc and more rapid cooling. Although 
such a feedback mode can be arranged to match the break in 
the present day luminosity function, a key issue is whether 
it explains the trajectory of star formation in galaxies over 
the past 5-10 Gyr. 

Similar progress has been made observationally in mea- 
suring the evolving stellar mass function of galaxies over 
0< z <1.5, where large and complete samples can be ob- 
tained (Fontana ct al. 2004; Drory, Bender & Hopp 2004; 
Bundy et al. 2006; Borch et al. 2006; Pozzetti et al. 2007). 
The advent of large-format near-infrared detectors used in 
conjunction with deep, spectroscopic and multi-wavelength 
surveys has characterized the evolving stellar mass function 
(Fontana et al. 2004), illuminated the bimodal nature of lo- 
cal galaxies (Kauffman et al. 2003), demonstrated the pres- 
ence of morphological evolution associated with assembly 
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since z ~ 1 (Brinchmann & Ellis 2000) and revealed how 
the quenching of star formation in massive galaxies produces 
the downsizing signature (Bundy et al. 2006). 

The time is therefore ripe for a direct confrontation be- 
tween recent simulations which incorporate "radio mode" 
feedback to fit the local luminosity function and the history 
of mass assembly over < z < 1.5 from the new generation 
of deep surveys. By using stellar mass estimates as a bridge 
between theory and observations we can gain significant in- 
sight not only into the success of the physical prescriptions 
employed by these models (Kitzbichler & White 2006), but 
also the utility of observations in answering the questions 
posed above. 

A key issue in comparing theoretical predictions with 
observations is the reliability of the latter. Previous com- 
parisons (Kitzbichler & White 2007; Bower et al. 2006; Cole 
et al. 2001) have been content to compare the total stel- 
lar mass function. However, the signature of downsizing in 
star formation is most readily tested by examining the mass 
function partitioned into star-forming and quiescent popu- 
lations (Bundy et al. 2006). While their fractional contribu- 
tions are relatively easily measured, the absolute numbers of 
these sub-populations are less certain because of the limita- 
tions of cosmic variance. This must be properly understood 
in any comparison, whether it be theory versus data or one 
dataset against another. 

As an example of the uncertainties arising from cos- 
mic variance, we note that while it is generally accepted 
in the community that the number of red-sequence galax- 
ies grows with time, confusion remains over whether this 
growth arises at the expense of a decline in the blue popula- 
tion (Bundy et al. 2006), or is dominated by "dry" mergers 
occurring within the red sequence, leaving the population of 
blue galaxies mostly invariant (Faber et al. 2007). In obser- 
vational samples broken by redshift bin and galaxy class, the 
impact of cosmic variance becomes a critical issue, sufficient 
perhaps to explain differences in empirical interpretations. 
In this paper we will use numerical simulations not only to 
test popular models of feedback against observations, but 
also to evaluate rigorously the limitations in the data aris- 
ing from cosmic variance. 

Our work follows logically from the earlier study of 
Kitzbichler & White (2007). Those authors considered the 
output of lightcones drawn from semi-analytical models in- 
corporating radio mode feedback as applied to the Millenium 
Simulation (Springel, White & Jenkins 2005) and compared 
predictions with various observables over < z < 5. They 
concluded broad agreement except in the abundance of high- 
z galaxies with large stellar masses. In this paper, we will 
focus on how star formation is distributed within the evolv- 
ing galaxy population. 

We construct an ensemble of 20 lightcones (j |3.3[ ) drawn 
from the Galform semi-analytic model (Bower et al. 2006) 
as applied to the Millennium Simulation. These lightcones 
are designed to mimic the Palomar/DEEP2 survey obser- 
vations presented in Bundy et al. (2006) which we argue 
in §2 below is currently the best dataset for addressing the 
questions of differential mass assembly and evolution of star 
forming and quiescent galaxies. These lightcones include the 
detailed survey geometry, optical and /sT-band magnitude 
limits, and photometric as well as stellar mass uncertainties 
described in Bundy et al. (2006). Critically, by studying the 



variation across our 20 realizations, we are also able to deter- 
mine reliable estimates of the effect of cosmic variance. As 
we will show, this is key to an accurate comparison between 
models and data. 

A plan of the paper follows. In §2 we justify our choice 
of the DEEP2/ Palomar spectroscopic stellar mass catalog 
(Bundy et al. 2006) as the comparison dataset. In §3 we 
discuss the Millenium Simulation and the associated Gal- 
form galaxy formation model. We review the various feed- 
back mechanisms and illustrate how they lead, in principle, 
to the concept of downsizing. We then describe how we con- 
struct suitable lightcones and implement the effect of mass 
errors in the context of the observational data. We compare 
the predicted colour distribution with the observations and 
discuss the uncertainties associated with various ways of se- 
lecting active and quiescent galaxies. In §4 we compare our 
mass functions with the observations of Bundy et al. (2006). 
We also discuss our findings with regard to how cosmic vari- 
ance may limit such conclusions, not only in the context 
of the current survey (Bundy et al. 2006) but also in other 
extant and projected surveys. 



2 OBSERVATIONAL DATA 

In this section we review the key features of the dataset 
presented in Bundy et al. (2006) which will serve as the 
"survey template" for constructing the lightcones described 
in §3 and the detailed comparisons discussed in §4. 

In order to derive reliable stellar mass functions, there 
are two highly desirable observational ingredients - exten- 
sive multi-band photometry extending to the near-infrared, 
and spectroscopic redshifts for the majority of the sample. 
Although no existing survey has complete spectroscopic cov- 
erage down to the faintest luminosities probed, we will focus 
in particular on comparisons of the downsizing signature, so 
good coverage at the high mass end is certainly advanta- 
geous. 

Our choice to focus on the DEEP2/Palomar catalog is 
motivated by a number of key advantages of this dataset. 
Foremost, as it builds on the extensive well-sampled DEEP2 
Galaxy Redshift Survey (Davis et al. 2003), it satisfies 
the two major criteria above. Furthermore, with respect 
to the important issue of cosmic variance, the survey cov- 
ers the largest area among published surveys with K-h&nd 
imaging — 1.5 deg^ spread over 4 independent fields. The 
C0MB017 survey (Wolf et al. 2003) covers ~0.8 square 
degrees in 3 fields, but has no near-infrared photometry. 
Although it has no spectroscopic coverage, it does benefit 
from highly-calibrated and well-tested photometric redshift 
data. The VVDS (Le Fevre et al. 2005) data set presented 
in Pozzetti et al. (2007) samples only a single field with an 
area of ~0.5 square degrees, and only ~0.2 square degrees 
has iC'-band imaging. 

We refer the reader to Bundy et al. (2006) for fur- 
ther details and only summarize the key aspects of the 
DEEP2/Palomar data set here. Palomar _K's-band photome- 
try was obtained in portions of all four fields targeted by the 
DEEP2 Galaxy Redshift Survey. The Extended Groth Strip 
(EGS) was the top-priority field and contains the deepest 
observations. In Bundy et al. (2006) three different redshift 
intervals were constructed, 0.4< z <0.7, 0.75< z <1.0, 1.0< 
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Table 1. A summary of the observational sample (Bundy ct al 2006). 



Field Number of Galaxies Area/deg^ 







0.4 < z < 0.7 


0.75 <z <1 


1 < z < 1.4 


Total 




EGS 




947 


1004 


722 


2673 


0.42 








209 


31 


240 


0.05 


DEEP2 < 






353 


266 


619 


0.13 








644 


409 


1053 


0.14 


Totals 




947 


2210 


1428 


4585 


0.74 



z <1.4 with T^s-band depths corresponding to Kab =21.8, 
22.0 and 22.2 for the three ascending redshift ranges. The 
lowest redshift sample comes entirely from the EGS, since 
galaxies with z <0.7 were excluded through colour selection 
in the other DEEP2 fields. In all cases, DEEP2 redshift tar- 
gets were limited to J? < 24.1. The field sizes are listed in 
Table 1. 

Central to our later discussion will be the partitioning of 
this sample into star-forming and quiescent galaxies. This is 
key to understanding the rate at which feedback suppresses 
star formation and provides valuable information in addition 
to the integrated stellar mass function (Kitzbichler & White 
2007). Bundy et al. (2006) considered both the rest- frame 
U ~ B colour and the rest-frame equivalent width of the [O 
II] emission line (Guzman et al. 1997) as proxies for the star 
formation rate. 



3 MODELLING 

Our aim in this work is to construct multiple sets of model 
galaxies selected in a similar manner, and with stellar masses 
inferred using similar techniques, as for the observed galax- 
ies. This will provide the fairest possible comparison between 
theory and data. 

3.1 Mass assembly 

The model galaxy samples are generated from the popu- 
lation of dark mater halos in the Millennium Simulation 
(Springel, White & Jenkins 2005). This simulation consists 
of approximately 10 billion dark matter particles each of 
mass 8.6 x 10*/i~^Mq evolving in a cubic volume of side 
500/i~^Mpc, assuming a ACDM cosmologjQ 

Dark matter halo merger trees are found from this 4- 
volume using the methods described by Harker et al. (2006). 
The lowest mass halos contained in these trees, of which 
there are about 20 million, consist of 20 particles corre- 
sponding to a total mass of 5 X lO'^ft-^Mo. Such halos could 
contain at most 9 x 10^h~^ Mq of baryonic material which 
is well below the lower limit of the stellar mass functions to 
be considered in this work. Therefore we do not expect the 
resolution of the Millennium Simulation to aflfect our results. 

The assembly of dark matter halos in ACDM is often 
described as "hierarchical" . This is appropriate in that some 
galaxies from one generation will merge to create the next. 

^ Specifically, a flat universe with = 0.045, Qm = 0.25, Hq = 
73km s-^Mpc-l and erg = 0.9. 
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Figure 1. The formation rate of galaxies within the Millennium 
Simulation as a fraction of their present day number. The up- 
per black line shows the total formation rate of galaxies with 
Mi, > IO^Mq , the lower blue line shows the rate at which pairs of 
these galaxies merge with each other. (Mergers involving smaller 
systems are not included.) The area between the curves will there- 
fore be equal to one. Redshifts relevant to the DEEP2/Palomar 
samples are indicated with vertical dotted lines. 

The importance of this contribution is illustrated in Figure 
[l] The detailed treatment of merging galaxies within a halo 
can be found in Bower et al. (2006). 

The Millenium Simulation demonstrates that the 
merger rate is expected to be quite small; less than 2% 
of galaxies with > 10® M© merge with each other ev- 
ery Gyr. At early times, minor mergers and the formation 
of new stars therefore create massive galaxies much faster 
than they be destroyed by mergers and the latter effect is 
almost negligible. As the universe evolves, the creation rate 
of new galaxies diminishes, leading to a near zero net growth 
in numbers (the difference between the two curves in Figure 
[l|. The growth trend is no longer in the number of galaxies 
but in their individual stellar masses. 

When one considers the differential mass assembly, 
therefore, mergers become particularly important. Figure [2] 
shows the total number of galaxies formed up to a given 
redshift through accretion (equivalent to the 'total' compo- 
nent in Figure [T]) alongside the number which remain after 
all mergers (since t=0) have taken place. The latter is the 
true, final number. 

The number density of galaxies in a given mass range 
will rise or fall depending on whether more galaxies arrive 
or leave that range due to an increase in their stellar masses. 
At intermediate masses, the population is almost unaffected 
by mergers; the creation and destruction rates are approxi- 
mately equal. At high mass, merging has a more significant 
effect because of the high ratio of potential progenitors to ex- 
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Figure 2. The evolving comoving number density of galaxies 
seen in the Millenium Simulation. Lines are coloured from black 
to red (top to bottom) corresponding to redshift limits {z=0.4, 
0.7, 1.0 and 1.4) adopted in the DEEP2/Palomar survey. The 
total number of systems formed by a given redshift is shown as 
a dotted line whereas the associated solid line shows the true 
number, including mergers. 



isting galaxies (indicated by a more steeply declining mass 
function at these masses). A particular deduction is that 
very massive galaxies (M* > 3 x 10^^ Mq) would be almost 
non-existent without merging. 

The preceding discussion is primarily based upon the 
formation and interaction of dark matter halos in the simu- 
lation. However, in order to construct Figure |2] we did have 
to consider the stellar content of the galaxies as well as their 
host halos. We now turn to discuss how stellar populations 
are introduced. 



3.2 Galaxies 

Halos in the merger trees introduced earlier are populated 
with galaxies using the Galform semi-analytic model, orig- 
inally described by Cole et al. (2000). Here we use the im- 
plementation of that model described in detail by Bower et 
al. (2006). The reader is referred to Cole et al. (2000) and 
Bower et al. (2006) for a full description. Here we briefly 
summarise the key physical processes included. 

The baryonic component of each dark matter halo is as- 
sumed to shock-heat during collapse of the halo to the virial 
temperature, at which point it settles into hydrostatic equi- 
librium and remains pressure supported until it can cool. 
The radius which has cooled at a given time after halo for- 
mation is calculated based on the metallicity of the gas, the 
cooling curve of Sutherland & Dopita (1993) and by assum- 
ing that it is isothermal at the virial temperature, with the 
radial density distribution: 



' core / 



(1) 



with rcorc set equal to one tenth of the virial radius. This 
choice is motivated by the simulations of Eke et al. (1998) 
and Navarro, Frenk, & White (1995). If the free-fall time 
for this radius has already been passed, any remaining en- 
closed halo gas is assumed to have settled into the central 
disk, where it becomes rotationally supported by its residual 
angular momentum and the combined gravitational poten- 
tial of the baryonic and dark matter. (The gas is otherwise 
added to the disk after a free-fall time) . 
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Figure 3. The fraction of galaxies with non-zero mass accretion 
due to cooling, plotted as a function of stellar mass, for four red- 
shifts. The bottom three span the range of the DEEP2/Palomar 
samples. Significant cooling with subsequent star formation oc- 
curs for all masses at high redshift, declining to be almost com- 
pletely suppressed for M« > IO^^Mq by z ~0.4. This demon- 
strates how the continued suppression of star formation due to 
radio mode AGN feedback can reconcile downsizing in the con- 
text of hierarchical models. 



We note that recent simulations have shown that, in low 
mass halos, shock heating may never occur, with gas arriv- 
ing into the central galaxy through "cold flows" (Keres et 
al. 2005; Birnboim & Dekel 2003). This mode of accretion is 
correctly accounted for in the model: In such low mass ha- 
los the cooling time will become much shorter than the halo 
dynamical time, and the mass infall rate will become inde- 
pendent of the assumption of shock heating and subsequent 
cooling and will instead be controlled by the cosmological 
mass accretion rate of the halo and the time required for 
free-fall. 

Star formation then proceeds in this cold disk gas with 
an instantaneous rate, ijj given by: 



^ ^* V200kms-^ 



Tdisk ' 



(2) 



where Tdisk refers to the dynamical time in the disk, Udisk 
the circular velocity at the half-mass radius and the dimen- 
sionless efficiency parameter, = 0.0029. 

The key necessity for feedback arises because the above 
assumptions lead to overly rapid cooling and star formation 
and hence to a galaxy mass distribution incompatible with 
observations (Benson et al. 2003). Various physical mech- 
anisms have been proposed to provide the necessary feed- 
back effects to prevent this from occurring. Two processes 
are particularly important: 

• Reheating of disk gas preferentially suppresses the for- 
mation of low-mass galaxies because of their shallow poten- 
tial wells. Some fraction of the disk gas is reheated by stel- 
lar winds and supernovae and is instantaneously (relative to 
the other timescales within the model) returned to the halo. 
This process is parameterised as follows: 



Mout = Pip, 



V '^disk / 



(3) 
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with Hhot ~ 485km/s. (3 can become extremely large, ex- 
ceeding 10* for small systems and approaching unity only 
for the most massive galaxie^ This effect cannot, therefore, 
be ignored when discussing the efficiency of star formation. 
In fact, the process of reheating as described by eqn. ([sjl is 
ultimately responsible for determining the eventual star for- 
mation rate and not, as might be initially expected, the star 
formation process as described by eqn. ([2|. This is discussed 
further in Appendix [A] 

• AGN heating was introduced by Bower et al. (2006) 
and occurs if the halo free fall time is shorter than the cool- 
ing time (formally, tcooi > dcooits, where Ofcooi = 0.58 is an 
adjustable parameter) such that a hydrostatic halo can exist. 
If this condition is met, further cooling of gas is prevented 
if the Eddington luminosity of the super massive black hole 
residing at the centre of the galaxy greatly exceeds the cool- 
ing luminosity (I/Edd > Esmbh-^'cooI, where esMBH = 0.04 
is an adjustable parameter). This strongly suppresses the 
formation of the most massive galaxies and imprints a near 
exponential cut-off in the abundance of the brightest galax- 
ies. The two free parameters, Ocooi and esMBH, were chosen 
by Bower et al. (2006) by constraining the model to match 
local luminosity function data. 

The black hole itself grows through gas accretion triggered 
by mergers and disk instabilities, acquiring Fbu = 0.5% of 
the available gas in each merger event and, similarly, 0.5% 
of the available gas in the galaxy if the disk's self-gravity 
exceeds the critical limit. 



GMdiek 
Tdisk 



> o.8v;^ 



(6) 



thereby causing the disk to become unstable. This criterion 
follows the work of Efstathiou, Lake & Negroponte (1982) 
though the particular constant, and the value of -Fbh, was 
found by requiring the model to match the Magorrian re- 
lation between bulge mass and black hole mass, Mbh ~ 
Mbuigo, as observed by Haring & Rix (2004). 

AGN heating is expected to primarily suppress the for- 
mation of the most massive galaxies. As a consequence of 
hierarchical growth, such systems will have considerable 
spheroidal components and, hence, massive central black 
holes. A study of AGN feedback implemented in semi- 
analytic models has been made by Croton et al. (2006), also 
using the Millennium Simulation, who found qualitatively 
similar results. 



^ These values for 13 axe extraordinary when viewed in the context 
of the physics thought to be involved in this outflow process. If 
eqn. (jsj is taken to be energy conserving, so that the gravitational 
potential energy gained by the expelled gas comes directly from 
the stellar winds and supernovae, the energy acquired by the gas 
per star formed would be 



3/8^ 



(485kms~^) 



(4) 



For a typical IMF approximately one supernova is expected for 
every IOOM0 of stars formed. This implies supernovae energies 
of about 



/200km/s 
V J^disk 



lO-*-*!. 



(5) 



This will exceed the available supply for galaxies smaller than the 
Milky Way, even without considering the efficiency with which 
this is conveyed to the gas (McKee & Ostriker 1977). 



Figure|3]illustrates the effect of AGN feedbaclj^in a simple 
way by showing the fraction of galaxies in the model with 
active cooling (i.e. those whose cooling has not been shut 
down by AGN heating). Cooling is largely unaffected in low- 
mass systems but is completely suppressed in the majority of 
massive systems. This illustrates how hierarchical structure 
formation can be made consistent with the observational 
phenomenon of 'downsizing'. 

In summary, in order for a galaxy to form a significant 
number of stars, there must be a supply of cooling gas from 
the surrounding halo to counteract the reheating of gas by 
the energy released. In this case equilibrium will quickly be 
reached which, under our particular parameterisation of the 
processes involved, results in a roughly constant specific star 
formation rat«[^ across all masses (see Fig. All. However, 
this balance is only possible in that fraction of galaxies for 
which the energy from the AGN is insufficient to prevent 
cooling of halo gas. In the case of no cooling, the fuel for 
star formation will quickly be exhausted, causing the galaxy 
to fade and redden. 

Thus, we can expect the galaxy population to be di- 
vided into two populations: those where cooling is occurring 
(high star formation rate, blue colours) and those with no 
cooling (little or no star formation, red colours). The relative 
abundance of each of these categories depends on the galax- 
ies' mass. The key question we wish to address is whether 
the growth and abundance of these two populations, which 
depends crucially on the feedback mechanism, matches that 
observed. 



3.3 Lightcones 

We use the prescriptions for gas cooling, star formation, 
and feedback described above to populate dark matter halos 
from the Millennium Simulation with model galaxies. The 
extremely large volume of this simulation allows us to reduce 
the statistical uncertainty in our predictions. However, to 
assess the significance of potential differences with observa- 
tions we must examine various sources of uncertainty. While 
errors based on photometry or Poisson statistics are typi- 
cally straightforward to estimate, the observational results 
of interest here — in particular, galaxy number density — are 
often dominated by sample or "cosmic" variance, which is 
much more challenging to quantify. The combination of a 
semi-analytic model and the Millennium Simulation pro- 
vides a powerful approach to this problem. After demon- 
strating general agreement between simulated and observed 
results, we construct numerous mock datasets with the iden- 
tical field geometry, magnitude limits, and photometric un- 
certainties as the real datasets. By comparing multiple re- 
alizations, we can accurately estimate the effect of cosmic 
variance. Below, we describe our method for constructing 
such lightcone dataset^ 



^ Figure 7 of Croton et al. (2006) illustrates a related point by 
showing the effect of AGN feedback in their model on the mean 
behaviour of all galaxies rather than the fraction which is affected. 

Specific star formation rate = M^/M*. 
^ Such datasets are referred to as "lightcones" (Kitzbichler & 
White 2007) as they contain all galaxies which intersect the past 
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Figure 4. Projection of a sample of galaxies on a lightcone drawn from the ^ lO^Mpc^ comoving volume of the Millennium Simulation. 
To mimic the DEEP2/Palomar survey, galaxies brighter than m^ AB = 24.1 are selected and split by star formation rate: star-forming 
(> O.2M0/yr; blue) and quiescent (< O.2M0/yr; red). The large black circle indicates the location of the observer. The lightcone is split 
into sections to fit on the page. Starting from the observer, the lightcones runs from left to right (with distance and redshift increasing 
along the lightcone). At the right hand edge of the page the lightcone reappears on the left hand side, shifted up the page by 60Mpc//i. 
Significant large scale structures — some of which span the width of the lightcone — are apparent, highlighting the need for a careful 
analysis of sample variance. 



We begin by populating the entire Millennium Simula- 
tion with galaxies using the methods described in §3.2[ This 
process is carried out for every available redshift in the Mil- 
lennium Simulation between z = 2 and z = 0. In addition to 
physical properties such as stellar mass and star formation 
rate we compute for each galaxy observable quantities such 
as apparent magnitudes in several bands and rest-frame ab- 
solute magnitudes. 

We then extract lightcones from this multi-redshift 
dataset using techniques similar to those of Kitzbichler & 
White (2007). The only significant difference lies in the 
method of interpolating galaxy properties between the out- 
put redshifts of the Millennium Simulation as discussed be- 
low. Briefly, we select random locations within the simula- 
tion volume and place an "observer" at that point. A line of 
sight is chosen following the methods of Kitzbichler & White 



lightcone of an observer located at some point within the simula- 
tion volume. 



(20070 and a cone with the geometry of the observed sam- 
ple is constructed around this line of sight. We then proceed 
to identify all galaxies which intersect the past lightcone of 
the observer within this cone. 

Kitzbichler & White compute magnitudes of each 
galaxy at its output redshift and the two adjacent output 
redshifts, thereby allowing them to interpolate to find the 
magnitude at any intermediate redshift. They effectively ap- 
ply a differential fc-correction to correct for the difference 
between the output redshift and the observed redshift at 
which the galaxy intersects the past lightcone. In our ap- 
proach, we track each galaxy from one output to the next. 



^ The lightcones constructed have an extent much greater than 
the size of the Millennium Simulation and so use the periodic 
boundary conditions of the simulation to effectively create a larger 
volume. The methods of Kitzbichler & White (2007) choose the 
line of sight in such a way as to minimize the possibility of the 
cone intersecting the same region of the simulation in different 
periodic replications. 
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Figure 5. The stellar masses of a representative 15,000 galaxies, 
derived from the simulated photometry, plotted as a fraction of 
the generated stellar mass. The IMF of Chabrier (2003) was used 
in the model and in the analysis used to make the mass estimates. 

Consequently, we can directly interpolate the properties of 
the galaxy to the precise redshift at which it is observed. Im- 
portantly, this allows us to include both fc-corrections and 
evolutionary corrections to magnitudes and to interpolate 
stellar masses which will change from one output to the next. 

Twenty lightcones are constructed, each having a solid 
angle of 0.42deg^, equivalent to the largest individual area 
covered by any of the three DEEP2/Palomar sub-samples. 
Galaxies are selected from within each lightcone if their -Rab 
and Ks magnitudes fall within the same limits as the ob- 
servational sample. The variation in number density across 
these 20 samples can then be calculated. 

We create additional sets of lightcones by trimming 
the originals to match the solid angles of the smaller 
DEEP2/Palomar fields which are relevant at 0.7 < 2: < 1.0 
(see Table [2|. For these redshifts, 4845 mock samples are 
generated, each consisting of four lightcones of the appro- 
priate size, and the variation measured across these sam- 
ples. The estimates of sample variance obtained in this way 
are therefore conservative; slight underestimates because of 
correlations between the different lightcones. 

3.4 Stellar mass estimates 

Stellar masses in the DEEP2/Palomar survey were de- 
rived for each survey galaxy using the spectroscopic redshift 
(when available) and a spectral energy distribution (SED) 
based on optical and near-infrared photometry. The effect 
of possible errors in this process were discussed in detail by 
Bundy et al. (2006). In the case of the simulated galaxies, we 
can evaluate the uncertainties independently by computing 
predicted photometric magnitudes using Galform and then 
applying the same method used by Bundy et al. (2006) on 
the observed photometry to rederive stellar mass estimates. 
Comparing these derived masses to the "input" values de- 
termined by the model provides a useful check on the stellar 
mass estimates and their uncertainties. 

For each of 15,000 galaxies, the predicted photometry 
for model galaxies was compared to actual data from the 
EGS field of the DEEP2 survey. Photometric errors were as- 



signed by randomly sampling the error distribution of EGS 
sources with similar magnitudes for each passband. In this 
way the perturbed magnitudes of the simulated galaxies re- 
flect the data quality of the EGS, including the variations in 
survey depth. The stellar mass was estimated as in Bundy 
et al. (2006) by comparing the "observable" SED of model 
galaxies to a large grid of stellar population templates and 
marginalizing over this grid, the final value being the median 
of the mass probability distribution. 

The comparison between the derived mass estimates 
and their input values are shown in Figure|5] For most galax- 
ies, the agreement is excellent with a scatter of ~0.15 dex 
as expected from the analysis in Bundy et al. (2006) of the 
internal uncertainties of the mass estimator and photomet- 
ric errors. This level of uncertainty is applied to the model 
masses when constructing mass functions below. Our mass 
comparison also identified a small subclass of model galax- 
ies where the rederived mass estimates are larger by ~0.1 
dex than their model values. These systems are of interme- 
diate mass and have star formation histories that have been 
sharply truncted by the Galform model, usually as a result 
of entering a larger halo. The timescales are typically much 
shorter than a dynamical time, leading to unphysical stellar 
populations that are poorly fit by the stellar mass estima- 
tor. Indeed, these galaxies populate a region of colour-colour 
space that is avoided by observed galaxies and represent the 
extreme end of the problem of over-quenching, which we will 
return to below. 



4 RESULTS 

Our primary goal is to compare the downsizing trends 
claimed by Bundy et al. (2006) via their colour-dependent 
stellar mass functions and to address the extent to which 
their evolving mass threshold is consistent with feedback 
models which reproduce the local mass and luminosity func- 
tion (Bower et al. 2006; Croton et al. 2006). 

4.1 Colour bimodality and the quenching of star 
formation 

We begin our comparison between the DEEP2/Palomar ob- 
servations and the Galform model with an analysis of the 
colour distribution. Figure |6] plots in three redshift intervals 
the restframe {U — B) colour-magnitude diagram of galax- 
ies drawn from one set of mock lightcones constructed to 
match the properties of the DEEP2/Palomar dataset. As 
discussed in Croton et al. (2006) and Bower et al. (2006), 
Figure [5] shows that the incorporation of radio-mode AGN 
feedback helps establish a distinct red sequence and pre- 
vents star formation in the brightest galaxies, as observed. 
The solid line in these diagrams has a slope given by van 
Dokkum et al. (2006) with a vertical offset chosen empiri- 
cally by Willmer et al. (2006) and Bundy et al. (2006) to 
divide the red sequence from the blue cloud in the observed 
DEEP2 distribution. This line has the equation: 

Ur - Br = -0.032 {Br + 21.5) + 0.454 - 0.25 (7) 

Clearly the colour cut defined by the observational data 
appears sub-optimal to split the model sample. A more ap- 
propriate cut (solid line) is therefore used in the investiga- 
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Figure 6. Galaxies from one of tlie 20 ligtitcone samples. As in 
Bundy et al. (2006), the star-forming and red-sequence galaxies 
are highliglited by applying a cut on the star formation rate at 
Af* = 0.2MQGyr~^. Galaxies above this cut are coloured blue 
while those below are coloured red in Figure |6] This criterion 
serves to illustrate the presence of two distinct populations but it 
does not identify each precisely. The solid diagonal line indicates 
the division made in the observational sample using rest-frame 
U — B colours. The solid line shows the division in the model 
sample, thereby producing Figures [s] & |9] 




Figure 7. The rest-frame U~ B colour distribution. Black points 
show observational values and shaded areas correspond to the 
total range of values found across 20 lightcones. For the latter, an 
observational error of 0.1 mag (1 s.d.) has been included. Both 
datasets are connected with smooth curves for visual clarity. The 
same star formation rate cut as Figurejo] M* = O.2M0Gyr~^, is 
applied to illustrate the division into two populations. 



tion into the evolution of these two populations (j ^4.2| and 
4.3 1. We note that this adjustement did not significantly 
change the relative fraction of red and blue galaxies, since 
the bulk of the red population remain above the dotted, 
observational line in Figure |6] 

A more profound discrepancy between the observed 
and modelled colours is shown in Figure [7| which plots the 
{U — B) distribution in the same three redshift intervals. 
Results from the lightcones are illustrated by shaded curves 
while data points signify the observed distributions from the 
DEEP2/Palomar dataset. Beginning with the total distribu- 
tion, we find that model galaxies trace a narrower range in 
{U — B) colour than observed galaxies. Even after includ- 
ing photometric uncertainties of 0.1 mag, the distribution of 
colours is simply too narrow, particularly at low redshift. 



Figure [7] reveals further insight into the problem that 
the Galform model has in reproducing the red sequence. 
The total distribution shows evidence for a red sequence in 
the two highest redshift bins, while for 0.4 < 2 < 0.7 the 
colour distribution looks nearly unimodal. Even in the two 
high-2 bins the red-sequence is offset blueward of the data — 
as we saw in Figure [6] — and seems to include larger numbers 
of red systems than observed. 
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Figure 8. The effect of cosmic variance on the derived stellar 
mass function at 0.4 < 2 < 0.7. Points are observational values 
from Bundy et al. (2006). Shaded areas show the 10-90% range 
of the simulated mass functions generated from 20 lightcones, 
defined as discussed in section f3.3| Black points and shading relate 
to the total mass function; red and blue to quiescent and star- 
forming components, divided as shown in Figure |6] Note that, 
in this figure, stellar masses are taken directly from the model, 
with no inclusion of errors arising from photometric uncertainties. 
The effect of including errors can be seen by comparison with the 
bottom panel of Figure [9] 



4.2 Evolution of the total mass function 

Despite the disagreement between tlie predicted and ob- 
served restframe colour distributions, we can still evaluate 
whether the current feedback prescriptions reproduce the 
general mass-dependent decline of star formation in galaxies 
since z ~ 1 which represents the basic 'downsizing' signa- 
ture. Such a signature should not be too dependent on the 
precise definition of what constitutes a quiescent or active 
galaxy. 

We begin by comparing the total mass function from 
both the mock lightcones and the observations. The results 
from the model are plotted in Figure [9] as the black shaded 
regions while the observations are indicated by black solid 
points. The width of the shading shows the 10-90% range of 
values recovered from the lightcone samples. As emphasized 
by Kitzbichler & White (2007), we find good agreement with 
the observed total mass function when the stellar mass un- 
certainties are convolved with the model results. This has the 
effect of increasing the high-mass end of the predicted mass 
function, as can be seen by comparing Figures [8] and [9] We 
note that we observe fewer low-mass galaxies than predicted 
at all redshifts, despite the extremely efficient conversion of 
supernova energy to ejected gas in this model ( ^3.2[ ). The 
completeness issues in the observational sample do not ex- 
plain this discrepancy because the observational magnitude 
limits have been included in the mock lightcones. 

Though disagreement between the predicted and ob- 
served total numbers appears to be minor, this thorough 
analysis of cosmic variance shows that there is still signifi- 
cant inconsistency with the data at many stellar mass inter- 
vals. 
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Figure 9. The evolving stellar mass function from the 
DEEP2/Palomar survey compared with the output of the semi- 
analytic models using the Galform radio-mode feedback. The 
predictions include the effect of errors in the stellar mass deter- 
mination. The 10-90% range of the various lightcone realizations 
are indicated by the shaded areas. All lightcones are reduced in 
size so that the four solid angles in each set correspond to those 
of the four DEEP2/Palomar survey fields (see Table [2]l. The sig- 
nificance of the colours is consistent with Figure |8] According to 
Bundy et al. (2006), the mass completeness limits from the Ks- 
band alone are log[M*/MQ] = 10.1,10.2 and 10.4 for the three 
respective redshift intervals. 



4.3 Verifying a decline in the star-forming 
population 

By adopting a suitable division in colour (Fig.[6|, the model 
galaxies are divided into red sequence and star-forming pop- 
ulations. It is clear from Figure [9]that the resulting compo- 
nent mass functions fail significantly to match the observed 
numbers, even after allowing for cosmic variance and mass 
errors. However, one important evolutionary trend is qual- 
itatively similar; the red-sequence population growing with 
time. The blue population in the model shows no significant 
evolution however, remaining static across the entire redshift 
range while the observed numbers fall with time. This leads 
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us to address an important discrepancy of interpretation in 
the observed mass functions. 

Bundy et al. (2006) argue that massive blue galaxies de- 
tected at z ~ 1 increasingly transform into quenched red sys- 
tems, leaving the total mass function relatively unchanged 
since z ^ 1. Faber et al. (2007) analyze B-band luminos- 
ity functions from DEEP2 and COMBO 17 and argue that 
the number of blue galaxies is essentially unchanging, with 
the increase in the massive red sequence population com- 
ing primarily from 'dry mergers' - mergers of similarly qui- 
escent galaxies. Bell et al. (2007) invoke a third option in 
which the red sequence is built primarily from quenched blue 
galaxies that are quickly replenished, leaving their numbers 
unchanged. The uncertain effects of cosmic variance have 
made distinguishing these scenarios challenging. With the 
improved measure of sample variance afforded by our mul- 
tiple lightcones, we can now evaluate the validity of these 
claims. 

The fundamental question is whether the evolution in 
the mass function of blue, star-forming galaxies claimed by 
Bundy et al. (2006) is significant given the uncertainties of 
cosmic variance. We have re-analyzed the observational re- 
sults of Bundy et al. (2006) to derive upper limits on the 
blue mass function at the high mass end for the two lower 
redshift intervals. These upper limits are determined by the 
value that would have been obtained had one blue galaxy 
been detected in each bin. In all cases, the upper limits are 
dn/dlog(M./M0) < lO-^'-^Mpc-^ 

As is apparent in Figure [9] this analysis reinforces the 
substantial decline observed in blue galaxies — nearly an or- 
der of magnitude across the redshift range at log(M,/M0) ~ 
11.4. This decrease is significant at the 99.9% level. The 
other mass bins above the _R-band determined completeness 
limit of log(Mt/M0) ~ 10.9 do not show significant evolu- 
tion, however. Using the upper limits, we find no evolution 
in the highest mass bin, log(M,/M0) — 11.7. It should be 
noted that no blue galaxies were detected at these masses 
with z < 1.0, so this result is a lower limit on potential 
evolution. At log(M,/M0) = 11.1, the observations show 
an increase from z~1.2toz~0.9 followed by an equally 
significant decrease from z « 0.9 to z ~ 0.5. Across the full 
redshift range, this is consistent with no evolution in the 
log(Af,/M0) = 11.1 mass bin at the 60% level. 

Taken together. Figure [9] and the analysis above sub- 
stantiates the claim made by Bundy et al. (2006) that the 
population blue, star-forming galaxies observed at 2: ~ 1 de- 
clines with time beyond the quenching mass, log(M,/M0) > 
11.2 (Bundy et al. 2006) . We note, however, that this evolu- 
tion is not easy to detect even in the large DEEP2/Palomar 
dataset. Although we see strong and significant evolution in 
one mass bin, averaging over the relevant mass range reduces 
the significance to the 2-3cr level after cosmic variance. 

Still, such evolution suggests that the red sequence is 
primarily built from the quenching of star-forming galaxies 
and only to a lesser extent from dry merging. A replenishing 
supply of blue systems would not seem to be necessary since 
such galaxies are not detected at lower redshifts. Finally, the 
substantial uncertainties from cosmic variance, which domi- 
nate the error budget, mean that much larger data sets will 
be required for more detailed studies capable of quantifying 
the importance of such effects as merging, internal growth 
due to SF, and transformations on the mass functions of 
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Figure 11. The colour magnitude distribution of all the galaxies 
in a subvolume of the simulation, generated using the model of De 
Lucia & Blaizot (2007). The solid black line indicates the colour 
divide applied to create Fig. |10[ 

different types of galaxies. Such work will represent an im- 
portant step forward in understanding the physical nature of 
the quenching mechanism. In Section [4.5[ these limitations 
are discussed in the context of upcoming surveys. 

4.4 Alternative models 

The method applied by Bower et al. (2006), populating the 
halos of the Millennium simulation using a galaxy formation 
model, was also adopted by Croton et al. (2006). The results 
from an updated version of this model (De Lucia & Blaizot 
2007) have been made publically avaiabl^ and the stellar 
masses and colours of all the galaxies in the simulation vol- 
ume were extracted. 

There are some significant differences in the way that 
some of the physical processes are modelled. Star formation 
is assumed to occur only if there is sufficient gas to exceed 
a critical surface density, following the work of Martin & 

^ http://www.g-vo.org/Millennium 
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Figure 10. Comparison between the stellar mass function predicted by the modelling techique and parameters of Bower et al. (2006) 
(dotted lines) and those of De Lucia &; Blaizot (2007) (solid lines). Both build on the halo population of the Millennium simulation. They 
are divided into two populations using the dividing criteria illustrated in Figures [6] and [TT] Points show the observational determinations 
of Bundy et al. (2006), but now with error-bars that reflect the cosmic variance calculated using the techniques of section [3. 3[ The points 
are from samples limited in magnitude at R=24.1 (except the highest redshift, which shows the photo-z supplemented mass function from 
Bundy et al. (2006), limited at R=25.1). Since the model mass functions are derived from instantaneous snapshots from the simulation 
(not from lightcones) the data is therefore included for reference, not for strict comparison. 



Kennicutt (2001). The mass converted into stars is then 5- 
15% of the gas which is above this value. This approach can 
be compared with equation (21. Gas is reheated from the 
disk in line with equation ([3 1 but with (3 = 3.5; the halo 
potential is involved only when calculating the gas ejected 
from the halo, not from the disk. 

AGN feedback is also incorporated, but rather than im- 
posing a cut-off black hole mass which prevents cooling com- 
pletely, the cooling rate is simply modified: 



M', = Mr, 



2 vir 



(8) 



The efficiency, rj — 0.1. Black hole accretion in this model is 



Mb 



relation mentioned in 



also set to match the Mbh 
section [3.21 

The stellar mass function derived from both models is 
shown in the left hand panel of Figure |10[ for the same red- 



shift ranges as before. It is worth emphasising that these 
number densities are not derived from lightcone samples, 
but from the entire volume. The points from the observa- 
tional sample are therefore included in Figure [To] for qual- 
itative companson only. Detailed quantitative analysis was 
possible using Figure[9] which sets the data against correctly 
constructed mock observations. 

Despite some considerable differences in their physical 
assumptions and choice of parameters, it is clear that both 
models produce very similar predictions. This similarity re- 
flects the fact that those parameters which remain theoret- 
ically uncertain are chosen to match each of the models to 
certain observational data. Similar data sets will have been 
used in both cases, and will have included the stellar mass 
function at 2: = 0. This partly accounts for the fact that the 
two models agree more closely with each other (and with 
the observations) at lower redshift. 
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In both cases, the predicted stellar mass function 
evolves in qualitative agreement with these observations. 
Major disagreement is restricted to distant, low mass galax- 
ies where survey completeness may provide an immediate 
explaination. 

The colour distribution of galaxies in this second model 



also shows a clear bimodality (Fig. Ill, as discussed in Cro- 
ton et al. (2006). The distribution is reproduced here to 
show the division between red and blue populations which 
has been used to produce the right hand panel of Figure 
|10[ showing their fractional abundance. As with the total 
number density, both models correctly predict the trend for 
the blue fraction to decreases with mass and increase with 
redshift, but specific values are frequently inconsistent with 
observations. 

The data points in Figure [lO] have updated error-bars 
which now reflect the cosmic variance estimates shown in 
Figure [9] The difference between the two models is rarely 
more than this cosmic variance, illustrating an inherent diffi- 
culty in constraining competing theories of galaxy evolution 
on the basis of their predicted stellar mass functions. 



4.5 The effect of cosmic variance on present & 
future surveys 

A salutory lesson from Figure[9]is that even with the largest 
and most complete survey of the evolving stellar mass func- 
tion to date, cosmic variance is still limiting quantitative 
conclusions concerning the growth rate of massive galaxies 
{M^, > 10" Mq) at the level of factors of order 2-3. These 
uncertainties have a larger impact than errors in the stellar 
mass (§3.4). 

Given we have developed the machinery to address this 
uncertainty, it is interesting to know the extent to which 
present and upcoming surveys will be limited by cosmic vari- 
ance in addressing basic questions: what is the assembly rate 
of the most massive galaxies from z < 2 to the present day? 

In an important paper, Somerville et al. (2004) intro- 
duced a simple analytical approach for estimating the effects 
of cosmic variance in a survey. These analytical estimates 
can be calculated for our mock samples by following the 
definition used in Somerville et al. (2004): 



2 

cry 



1 

(TV)- 



(9) 



It can also be calculated using the correlation length, which 
is available for the DEEP2 survey (Coil et al. 2008). Fig- 
ure [12] compares this analytical prediction with the variance 
derived directly from our multiple lightcones technique. As 
can be seen, the relative agreement is quite encouraging. 
Returning to the physical questions posed above. Figure 
shows the number density of intermediate-mass (10^*^ 

11\ 1 u:„l, ™„„„ ^ inll 



13 



< 

M/Mq < 10'^) and high -mass (M > IO^Mq) galaxies for 
the entire volume of the Millennium Simulation. The key 
question is the extent to which present and upcoming sur- 
veys can verify not only the basic trends of mass assembly 
but also the mass-dependent variations. 

First, we consider the situation for the data evaluated in 
the present paper. For the Bundy et al. (2006) data plotted 
in Figure [l3l errors correspond to the 10-90% range of values 
found from the 20 lightcone realizations of this survey. Over 
0< z <1.4, the predicted growth trends in intermediate mass 
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Figure 12. The relative cosmic variance of the red and blue 
galaxy populations as a function of survey volume. The solid lines 
are calculated using the 20 lightcones, the dotted lines using the 
analytical approach of Somerville et al. (2004) (see text). The 
relative variance of the total population of galaxies across the 20 
lightcones is shown in black. 



Table 2. Some details of the selection of observational surveys 
referenced in Figure [13] 



Field 




Arca/deg^ 


Mag 


Limits 








Rab 




EGS & DEEP2 




Table 1 


24.1 


20.5 


GOODS North 




1 X 0.042 


24.1 


21 


COSMOS 




1 X 1.5 


25 


22 


UKIDDS ^ 


D.E.S. 
U.D.S. 


4 X 8.75 
1 X 0.77 


23 
27 


21 
23 



galaxies (upper, black line) are only marginally confirmed 
and the predicted growth in massive galaxies (lower, red 
line) is inconsistent with observations. 

Much of the differential trends observed are affected by 
the area of the survey and the optical magnitude limit (Ta- 
ble 2). The situation with earlier HST samples (GOODS, 
Giavalisco et al. (2004)) is considerably worse. However, 
even with upcoming panoramic surveys underway with HST 
(COSMOS, ScoviUe et al. (2007) ) and UKIRT (UKIDSS, 
Lawrence et al. (2007)), the situation is not significantly bet- 
ter. Thus precise observational constraints concerning the 
differential evolution of mass assembly must await a future 
generation of surveys. 
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Figure 13. The predicted comoving number density of galaxies derived from the full volume of the Millenium Simulation for intermediate- 
mass galaxies (black curves) and those with high stellar mass (red curves). Shaded areas indicate the constraints on mass assembly in 
these two mass bins according to the results of Bundy et al. (2006). The mean number densities for 10< log(M*/MQ) <11 may be affected 
by some incompleteness, which becomes important at log(M*/MQ) < 10.4 in the photo-z supplemented sample used here (Rab < 25.1). 
However, the mock observations confirm that most of the discrepancy between the data and the models at lower masses is real and not 
simply the result of incompleteness (Fig|9]l. The shading indicates 10-90 % error due to cosmic variance, calculated using the techniques 
described in i ]3.3| Error bars illustrate the equivalent error expected in other present and upcoming surveys (left to right are listed from 
top to bottom). For convenient comparison, these errors are all given relative to the same, observationally determined mean densities. 
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5 SUMMARY 

Lightcones derived from the Millenium Simulation, designed 
specifically to match the geometry and selection parameters 
of the DEEP2/Palomar survey of Bundy et al (2006), re- 
produce reasonably well the evolving stellar mass function 
over 0.4< z <1.2. These mock observations are populated 
with galaxies using the Galform code, which incorporates 
a prescription for 'radio mode' feedback. Such feedback pro- 
cesses are thus an adequate explanation of the broad trend 
of 'downsizing' seen in recent redshift surveys. 

However, the Galform code can not satisfactorily re- 
produce the evolving form of the bimodal colour distri- 
bution, or the stellar mass functions of the corresponding 
two populations. An alternative model (De Lucia & Blaizot 
2007), proposed originally by Croton et al. (2006), produces 
similar discrepancies; both models appear to over-quench 
star formation in intermediate mass galaxies. 

This approach also allows us to address the important 
question of how cosmic variance may limit the validity of 
various conclusions drawn from observational surveys. Us- 
ing variance estimates derived from the lightcones, we con- 
firm the significance of the decline since z ~1 in the num- 
ber density of massive blue galaxies claimed by Bundy et 
al (2006). We argue that the transformation of these blue 
galaxies must, necessarily, provide the bulk of the associ- 
ated growth in red sequence galaxies given the near-constant 
total mass function. 

This discussion of cosmic variance is extended to 
demonstrate the limitations of other, more ambitious, sur- 
veys ongoing at the present time in terms of detecting the 
mass-dependent growth of galaxies since z ~1. 
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APPENDIX A: STAR FORMATION 
EFFICIENCY 

In this Appendix we illustrate the importance of gas ejection 
from the disk, which was given in eqn. ([3|, for determin- 
ing the steUar masses of galaxies. The differential equation 
which controls the mass of cold disk gas, Mcoid, is 



Afc, 



(Al) 



where T^, = rdisk/e*/(i'dis5k/200km/s)^/'^. The parameter 
R = 0.39 is the fraction of material in stars which will be 
recycled back into the interstellar medium in the instanta- 
neous recycling approximatiorj^ 

For the relevant range of circular velocities, 5 < /3 < 
10, 000, so outflow will be fierce and the only way that a 
significant supply of cold gas can be retained in the disk is 
through constant replenishment from the cooling halo gas. 
This corresponds to the limit A/coid for which equation 
(|ATJ has the following solution: 



Tcft 



l-R + 13 



(A2) 



Given the large values of (3, the effective timescale, t^s will 
be extremely short, leading quickly to a steady state where 
the rates of cooling and reheating gas are approximately 
constant and equal. In this case, cold gas mass and star 
formation rate are given by: 



Mcold ~ MinTeff, M^^Min- ^r--—. (A3) 

1 — H -T p 

(Note that the star formation rate given here is the rate of 
increase of mass in long-lived stars, hence the inclusion of 
the factor 1 — R.) Because of the tendancy towards these 
limits, the star formation rate in any particular galaxy is 
approximately constant (assuming a constant Min), ignoring 
short term perturbations resulting from mergers and grav- 
itational instabilities which contribute little over cosmolog- 
ical timescales. In this approximation, the stellar mass is 
given by the star formation rate integrated over the avail- 
able time, t{z). 

The specific star formation rates for galaxies at redshift 
z should therefore be expected to be scattered about the 
value 



M^{z) 



1 



(A4) 



Figure |A1| demonstrates that these analytic approxima- 
tions are broadly consistent with the calculated values for 
~ 400, 000 galaxies in a sub- volume of the model. 

The above argument implies that specific star forma- 
tion rate is independent of the equations which explicitly 
govern it. The parameters for star formation efficiency, in 
eqn. ([2|, and the strength of outflow, /3 in eqn. (|3| will have 
only a secondary influence of the underlying trend described 
by eqn. ((A4|). 



° This particular value is calculated from the initial mass function 
of Kennicutt (1983): 
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Figure Al. An iUustration of the relationship between cooling 
rate and star formation rate as a function of disk circular speed, 
Va for ^ 400, 000 galaxies in an arbitrary sub- volume of the model 
at redshift z = 0.5. The top panel shows the (weighted) fraction 
of galaxies which have any cooling at all and so appear on the 
other three sets of axes. The second panel shows the cooling rate 
itself and the third panel shows both cooling rate and the star 
formation rate as a fraction of the existing steUar mass of each 
galaxy. The dashed line shows the constant star formation rate 
which would be required to accumulate the relevant stellar mass. 
The bottom panel illustrates that the ratio between the star for- 
mation and cooling is a function of the outflow efficiency, f3{vc), 
as predicted by eqn. (|A3[|. 



The physical property which is controlled by the star 
formation efficiency is the mass of cold gas in the disk rel- 
ative to the supply from the halo and, consequently, the 
stellar mass: 



t{z) ■ 



(A5) 



So, more efficient star formation leads to a lower gas fraction 
but not to higher star formation rates for galaxies of a given 
mass. 

The cooling rate and outflow efficiency, 13, determine 
the size of the halo in which a galaxy of a certain mass 
must reside in order to fuel its star formation. This is how 
the smaller-scale physics of star formation and supernovae 
infiuence the collective properties of galaxies on cosmological 
scales. 
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